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SURVEY ON THE INVERSE SPECTRAL PROBLEM 


STEVE ZELDITCH 


Let (M, g) be a Riemannian manifold of dimension n, possibly with boundary, and denote 
its Laplacian by 


d Xj : 


where g t j = g(-£g, gf-), [ g y ] is the inverse matrix to [g ig ] and g = det [g l3 \. When M is 
compact the spectrum of A g is the discrete set of eigenvalues 


-Ag(pj = A foj, (<Pi, ifj) = Sij 

Bipj = 0 on dM, 


( 1 ) 


repeated according to their multiplicities. When M has a non-empty boundary dM, one 
imposes boundary conditions such as Dirichlet Bu = u\qm or Neumann Bu = 8 v u\qm- 

More generally, we may consider the eigenvalue problem for semi-classical Schrodinger 
operators 

Hh'dh.j ^ 2" A + V’fij = (2) 

where V G is a potential. Here, h is the Planck constant. When ( M,g) is non¬ 

compact, the spectrum of A g is often continuous, and one might instead pose the inverse 
problem for eigenvalues of the scattering operator Sn (‘phase shifts’) or poles of its ana¬ 
lytic continuation (‘resonances’) 1J6J If the potential V grows at infinity, H\ has a discrete 
spectrum. 

The inverse spectral or scattering problem is to determine as much as possible of the metric, 
potential or domain from spectral or scattering data. In the case of discrete spectrum one 
would like to “invert” the map Sp : (. M,g,V ) —>• Spec (Hn) = {Aj}°‘L 1 , or at least to show 
that it is 1-1 when restricted to a reasonable class of ( M,g , V). We note that there are two 
types of inverse problems for Schrodinger operators ([2]): (i) the more difficult one where h is 
fixed, and (ii) the substantially simpler one where one attempts to determine the potential 
from the 1-parameter family of spectra {Ej(h)} < ^ =1 . In some sense, one would like to use 
the A j as “coordinates” on the ‘space’ of isometry classes of metrics, potentials or domains. 
This is only a heuristic way to envision the problem because of its infinite dimensionality. 
Except in dimension one, the range of Sp is not known and is certainly a very small subset 
of the space £ C of sequences compatible with the Weyl law for counting eigenvalues, 
which in the purely metric case takes the form #{j : A j < A] ~ C n Vol(M, g) X n . However 
it leads one to ask if Sp is ‘generically’ 1-1, to study its first and second derivative, to 
describe its level sets (isospectral sets), or to see if there is any natural additional data one 
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may use to supplement Sp to make it a well-behaved map. A simpler problem is spectral 
rigidity: can one deform a metric, potential or domain in a non-trivial way while keeping 
the eigenvalues fixed? I.e. do there exist (smooth) curves in an isospectral set of domains, 
metrics or potentials with a fixed spectrum? Isometric domains or metrics are regarded as 
the same, or as trivially isospectral, and when we discuss isospectral domains or metrics it 
is understood that they are non-trivially isospectral. 


0.1. Historical background. The inverse spectral problem was first stated by Sir Arthur 
Schuster in 1882 (see Appendix F of [St] (p. 395)) :“...Find out the shape of the bell by means 
of the sound which it is capable of sending out. And this is the problem which ultimately 
spectroscopy hopes to solve in the case of light. In the meantime we must welcome with 
delight even the smallest step in the right direction.” An early result is that of Rayleigh, 
determining the coefficients of a differential operator from its eigenvalues [R]. The problem 
is best known in the formulation of L. Bers, S. Bochner, and M. Kac (“Can you hear the 
shape of a drum”). Here, the metric is assumed to be Euclidean, V = 0 and one wants to 
determine a domain C M 2 from its Dirichlet or Neumann eigenvalues. 

In quantum physics, the inverse problem is mainly to determine the potential (or metric) of 
(Ej) from it spectral and scattering data |Nl !N2i IChS] , The idea is that this data is observable 
in experiments while the potential is not. There are many types of scattering data but in this 
survey we only discuss the phase shifts or eigenvalues of the scattering matrix Sh (a unitary 
operator acting on the sphere; dH]) • The problem of recovering the potential of a Schrodinger 
operator from eigenvalues and phase shifts was studied since the 1940’s when V. Bargmann 
discovered two potentials with the same phase shifts [Bari , Bar2j. The inverse problem in 
dimension one was solved with some restrictions by Gel’fand-Levitan and Marchenko (see 
[GLllM arl lMar2~] and also |Dy[ |Dy2| for background and history). Some other early articles 
are eubiIv n. 

The inverse spectral problem in one dimension is simpler than and quite different from 
the problem in higher dimensions and in this survey we only consider the problem when 
dim M > 2. To give an idea of how the difficulty grows with the dimension, it is simple to 
prove that the standard metric g 0 on the sphere S n is determined by the eigenvalues of A go 
when n = 2. It is difficult but possible to show g 0 is determined by Spec(A go ) for n < 6 
(Tanno). It is unknown whether g 0 is determined by Spec(A go ) when n > 7, although they 
are known to be spectrally rigid and locally determined (Tanno). On the other hand, it is 
known that balls in M n are determined by their eigenvalues (by the isoperimetric inequality). 

There exist two broad classes of results in inverse spectral theory: (a) construction of 
isospectral pairs or deformations, i.e. of ‘counter-examples’ to the inverse spectral problem; 
(b) positive results proving rigidity or spectral determination on classes of domains, metrics 
or potentials. In this survey we concentrate on positive results of type (b) and ask, how 
much of the geometry is determined by the spectrum? We refer to ea for some results of 
type (a). 

There are a number of recent surveys on the inverse spectral problem [ DH!, Zlj . The 
present survey is shorter and less systematic. Its purpose is to give a reader who is not 
very familiar with inverse spectral theory an idea of the current state of the art as well as 
some idea of what is involved in proving inverse spectral results. Some are new and some 
are relatively old but illustrate different features of the problem. In particular, we do not 
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attempt to provide a comprehensive list of references (see [DHUzU). One of the charms 
of the subject is that the problems have an obvious motivation in both mathematics and 
physics, are easy to state, yet remain largely open after 50 years. After reviewing the basic 
techniques we state some tantalizingly simple sounding problems in 

0.2. Analyticity, symmetry, accidental degeneracies and multiple spectra. The 

inverse spectral problem in higher dimensions involves a vast collection of objects and it is 
not surprising that there are few general results applying to all of them. Inevitably there 
exists sporadic spectral behavior that allows isospectral sets to have surprising properties. 
For instance, the eigenvalues or lengths of closed geodesics might have unusual multiplicities. 
In the theory of automorphic forms, the Langlands (-Eichler) correspondence between arith¬ 
metic hyperbolic surfaces with cusps and compact arithmetic surfaces defined by quaternion 
algebras is a (partial) isospectrality phenomenon that has no known global geometric inter¬ 
pretation. The correspondence is proved by ‘matching orbital integrals’, e.g. by comparing 
lengths of closed geodesics on both surfaces. Thus, the partial isospectrality has arithmetic 
rather than geometric origins. 

It is common to add assumptions or restrictions in order to prove positive results. Before 
introducing definitions and results, we call the reader’s attention to some common ones. 

• (A) The potentials, metrics or domains are often assumed to have a symmetry, either 
discrete or continuous. On M ra , potentials or domains are assumed to be invariant 
under reflections Xj —Xj , or under a di-hedral symmetry |Z2j . Continuous symme¬ 
try often takes the form of radial potentials, surfaces of revolution etc. The inverse 
spectral problem is restricted to the class of such symmetric domains or potentials. 

• (B) The metrics, potentiasl or domains are often assumed to be real analytic. Convex¬ 
ity is another natural condition. In the counter-example of [GWWj of non-isospectral 
domains with the same Dirichlet (resp. Neumann) spectra, the domains have corners 
and are non-convex. 

• (C)The isospectrality is often assumed to occur for more than one spectrum, as in 
the case of ([2]) for the one-parameter family h —>- Ej(h ) of spectra. Or in the case of 
domains, one might prescribe both the Dirichlet and Neumann spectra, or both the 
spectrum ‘inside’ the domain and also ‘outside’ of it. 

• (D) In counter-examples, the metric, potential or domain often has some non-generic 
behavior such as high-dimensional manifolds of closed geodesics or other mulitplicities 
in the length spectrum. Most positive results require multiplicity one in the length 
spectrum. 

Let us explain the background to these assumptions. 

In the early days of one-dimensional inverse spectral theory, Borg [Boj, Levinson [L] and 
Marchenko |Mari IMar2) proved that a potential V G L 1 [0, L] is determined by two spectra 
of a Schrodinger operator (somewhat similar to prescribing spectra of ([2]) for two values of 
h). When the potential is symmetric, q(L — x) = q(x) then one good spectrum, either the 
Dirichlet or Neumann spectrum, suffices (see e.g. [BrH] ). The solution is not unique for 
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general non-symmetric potentials. Results assuming Z 2 symmetries are multi-dimensional 
generalizations of this result. See |Z2l IGU, iHel IHZ2I] for recent results of this kind. 

Radial symmetry is a very common assumption in inverse scattering theory since it es¬ 
sentially reduces the inverse problem to one dimension [KKS> iN3] . The author is not aware 
of any inverse results for phase shifts that do not assume radial symmetry. The analogous 
problem for metrics is to assume the surface has an S 1 rotational symmetry, i.e. is a surface 
of revolution. When it is also mirror symmetric, the inverse problem was solved in [ BrH] by 
reducing the inverse problem to that of Borg-Levinson-Marchenko. There are several other 
results of this kind due to M. Kac, D. Stroock and P. Berard (see [DHJ). In |Z3] the author 
solved the inverse spectral problem for convex analytic surfaces of revolution without mirror 
symmetry; the proof used what we now call a normal forms approach (Theorem 13.41 and 

S3}- 

A plane domain may be locally defined as the graph of a function y = f(x) over (—1,1) and 
therefore has the same number of functional degrees of freedom as a surface of revolution or 
a radial potential. I.e. the unknown is a function of one real variable. If the domain has the 
Z 2 x Z 2 symmetries of an ellipse, the problem has the same symmetries as the one dimensional 
Schrodinger problem of Borg-Levinson. The lack of uniqueness in the 1 D case might suggest 
lack of uniqueness for the plane domain problem without two symmetries. However, if the 
domains are real analytic, the one Z 2 symmetry is enough to solve the inverse problem for 
the Dirichlet or Neumann spectrum ([Theorem 13.11 |Z2j ). It is not known at present if the 
inverse spectral problem is uniquely solvable (even in the analytic class) without assuming 
at least one symmetry. 

The history of results for radial potentials or surfaces of revolution might also suggest that 
one can only solve the inverse problem if the unknown is a function of just one real variable, 
e.g. that one cannot ‘hear the shape of a multi-dimensional drum’. However this is also 
not true. I 11 [HZ ] H. Hezari and the author proved that one can determine a real analytic 
ZJJ-synnnetric domain in M n for any n from its Dirichlet or Neumann spectrum (Theorem 
13.3p : here the unknown is a function f(x\,... ,x n ) of n variables. An earlier result for the h 
family of spectra of (J2J) for Z 2 -symmetric potentials was proved in [GUI [He]; in [ He] fewer 
symmetry assumptions were needed. 

Real analyticity is often assumed because the spectral data often leads to local expressions 
in terms of Taylor coefficients of the metric, potential or domain around closed geodesics. 
Unless one is making a deformation, there is no obvious way to relate spectral information 
for distinct geodesics. Real analyticity is used to make a local result global, i.e. to use Taylor 
coefficients at one point to determine the unknown function. 

The //-problem is significantly simpler than the problem for a single spectrum. It is similar 
to assuming the knowledge of the spectrum of several commuting operators, e.g. the joint 
spectrum of A g and of L — on a surface of revolution; here L generates rotations around 
the axis. I 11 [Z2], the principal step is to show that the single spectrum of A g determines the 
joint spectrum. This is clearly necessary if the inverse spectral problem for one value of h is 
solvable. 

Regarding (D), we recall that isospectrality of two Laplacians A j is the condition that 
there exists a unitary operator U such that UA\U* = A 2 . In the correspondence princi¬ 
ple between classical and quantum mechanics, the classical condition is that there should 
exist a symplectic diffeomorphism % between the cotangent bundles which restricts to a 
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diffeomorphism x '■ S gi M\ —* S* 2 M 2 between the unit cosphere bundles and which conju¬ 
gates the geodesic flows G f - in the sense that \G\x~ l = G\. The correspondence principle 
only works if U is a special type of operator known as a Fourier integral operator. We 
refer to | Zw ] for background. In [ Z5] it was pointed out that the symplectic maps underly¬ 
ing unitary Fourier integral operators can be multi-valued correspondences, and that all of 
the well-known counter-examples of Sunada [SuJ were conjugate by unitary Fourier integral 
operators (the unitary conjugation was also observed by Berard, who termed them trans¬ 
plantations). Thus, ‘most’ of the counterexamples are Fourier-isospectral but there are some 
examples which are not. 

As discussed in 1 14.31 one of the main approaches to inverse results is through Birkhoff 
canonical forms, which are special ‘local’ unitary Fourier integral conjugations between 
Laplacians in small neighborhoods of closed geodesics. This is a different kind of condi¬ 
tion from local isometry of the metrics but is related to it in special cases. See |Gor \m 
for discussion of isospectral pairs which are not locally isometric. 

As this discussion indicates, a natural direction for future work in inverse spectral theory 
is to remove the symmetry and analyticity assumptions as much as possible. It would be 
useful to know if (or to what extent) symmetry or analyticity were a spectral invariant for 
some reasonable classes of metrics, potentials or domains. The only result of this kind known 
to the author is the recent one of II)1 IV . which shows that, in the case of the 1-parameter 
h inverse spectral problem, a radial monotonic potential Vo of a Schrodinger operator (J2]) on 
IR n is spectrally determined among ‘all’ smooth potentials (not assuming any symmetry of 
the latter). 

With this background in mind we now present the basic objects in inverse spectral theory. 


1. Spectral invariants 

The data from which one hopes to recover the metric or potential is the list {A^}“L 0 of 
eigenvalues, enumerated with multiplicity, i.e. the dimension of the eigenspace ([I]). We first 
consider how to assemble the data into useful invariants. 

1.1. Traces and trace formulae. Since the origins of inverse spectral theory, the principal 
spectral invariants are defined by trace formulae. One assembles the eigenvalue data into 
generating functions such as the heat trace 0(f) := e ~ tXj or wave trace 

S(i) :=]T e ia J. (3) 

3 

The simplest way to relate spectrum and geometry is to study the singularities of these 
spectral functions. 

The heat trace is smooth except at f = 0, and the coefficients of its expansion at f = 0 
were the first spectral invariants to be studied. In dimension 2, the heat trace has the form 
0(f) = f~ 1 Area(M, g ) + ^y(M) + A f M r g dV g H— •. There is just one singular term and its 
coefficient is the integral of a local geometric invariant. One also sees that all of the Taylor 
coefficients of fO(f) are given by integrals of local geometric invariants. This shows that 
there are useful spectral invariants besides the ones which arise from singularities. 
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1.2. Length spectrum and wave trace singularities. The length spectrum of a bound¬ 
aryless manifold (M, g ) is the set 

Lsp(M, g) = {L 71 < L 72 < • • • } (4) 

of lengths of closed geodesics 7 j, i.e. the set of distinct lengths, not including multiplicities. 
One refers to the the length spectrum repeated according to multiplicity as the extended 
length spectrum. 

For a generic metric on a manifold without boundary, the length spectrum is a discrete set 
and moreover the lengths have multiplicity one. In the boundary case, the length spectrum 
Lsp(Q) is the set of lengths of closed billiard trajectories, which has accumulation points at 
lengths of closed billiard trajectories which glide for some interval of time along the boundary. 
In the case of convex plane domains, e.g., the length spectrum is the union of the lengths of 
periodic reflecting rays and multiples of |<9f2|. 

1.3. Singular support of the wave trace. The wave group of ( M,g ) is the unitary group 

U(t) = : L 2 (M) ->• L 2 (M). (5) 

Its trace ([3]) is a distribution on M. The first result on the wave trace is the Poisson relation 
on a manifold without boundary, stating that the singular support of S(t) is contained in 
the length spectrum, 

Sing Supp S(t) C Lsp(M,g). ( 6 ) 

It was proved by Y. Colin de Verdiere CdV3] , Chazarain |Ch2j . and Duistermaat-Guillemin 
[DG ] (following non-rigorous work of Balian-Bloch [ BB2] and Gutzwiller). The generalization 
to manifolds with boundary was proved by Anderson-Melrose [ AM] and Guillemin-Melrose 
[GMj. As above, we denote the length of a closed geodesic 7 by L 1 . For each L = L 1 G 
Lsp(M,g ) there are at least two closed geodesics of that length, namely 7 and 7 _1 (its time 
reversal). The singularities due to these lengths are identical so one often considers the even 
part of S(t) i.e. TrE(t) where 

E(t) = cos(f a/—A 9 ). (7) 

The containment relation (J 6 ]) could be strict if a length L G Lsp(M, g) is multiple. In this 
case, S(t) might be smooth at L G Lsp(M,g). 

One cannot determine multiplicities of lengths from the singularities of S(t) in any simple 
way. The Sunada-type isospectral pairs |Su] always have multiple length spectra, and for 
many (presumably, generic) examples, the length spectra have different multiplicities. 

In the notation for Lsp(M, g) we wrote L 7j . as if the closed geodesics of this length were 
isolated. But in many examples (e.g. spheres or flat tori), the geodesics come in families, 
and the associated length T is the common length of closed geodesics in the family. In place 
of closed geodesics, one has components of the fixed point sets of G T at this time. 

1.4. Singularity expansions. The next fact is that S(t) has a singularity expansion at 
each L G Lsp(M,g ): 

S(t) = e 0 (t) + Y,LGLsp(M,g) m ° d C °°, (8) 

where eo, are Lagrangean distributions with singularities at just one point, i.e. sing supp eo 
{0}, sing supp eL = {L}. When the length functional on the loopspace of M is a Bott-Morse 
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functional, the terms have complete asymptotic expansions. In the Morse case (i.e. bumpy 
metrics), the expansions take the form 


e 0 (t) — ao t - n (t + iO) n + Uo,-n+i(^ + *0) " +1 + ■ • • 


(9) 


eL(t) = a L -i(t — L + iO) 1 + a L ,o log(t - (L + i 0 )) 

( 10 ) 

+ a Lj i(t - L + iO) log(t - (L + i 0 )) H - , 

where refers to homogeneous terms of ever higher integral degrees ([DG]). The wave 
coefficients a at t — 0 are essentially the same as the singular heat coefficients, hence are 
given by integrals over M of f M Pj(R,VR, ...)dvol of homogeneous curvature polynomials. 
The wave invariants for t ^ 0 have the form: 


a L,j 


, 3 1 

7 \Lry=L 


( 11 ) 


where a 7J - involves on the germ of the metric along 7. Here, {7} runs over the set of closed 
geodesics, and where L 7 , L#, m 7 , resp. P 7 are the length, primitive length, Maslov index and 
linear Poincare map of 7. The Poincare map is defined in ((131) . For instance, the principal 
wave invariant at t = L in the case of a non-degenerate closed geodesic is given by 


dL -1 


£ 


e'^-'L* 

det(/-P 7 )|5 


( 12 ) 


The same formula for the leading singularity is valid for periodic reflecting rays of compact 
smooth Riemannian domains with boundary and with Neumann boundary conditions, while 
in the Dirichlet case the numerator must be multiplied by (—l) r where r is the number of 
reflection points (see [ GMl I'PS] ). 

The wave invariants for t 7 ^ 0 are both less global and more global than the heat invariants. 
First, they are more global in that they are not integrals of local invariants, but involve the 
semi-global first return map of the closed geodesic 7T,. One could imagine different local 
geometries producing the same first return map. Second, they are less global because they 
are determined by the germ of the metric at 7 and are unchanged if the metric is changed 
outside 7. 

Thus, associated to any closed geodesic 7 of ( M,g ) is the sequence {a 7 r j} of wave invari¬ 
ants of 7 and of its iterates Y. These invariants depend only on the germ of the metric at 
7. The principal question of this survey may be stated as follows: 


Problem 1 . How much of the local geometry of the metric g at 7 is contained in the 
wave invariants Can the germ of the metric g at 7 be determined from the wave 

invariants? At least, can the symplectic equivalence class of its germ be determined? 

As will be discussed below, the classical Birkhoff normal form of the metric (or the Poincare 
map P 7 ) at 7 is determined by the wave trace invariants [Gj (G2( |Z3j (Z4j . 

1.5. Domains and metrics with the same wave invariants. Although we are empha¬ 
sizing wave invariants, it is known that they do not generally determine a domain or metric. 
There exists an example known as a Penrose mushroom (due to Michael Lifshits) which 
shows that wave invariants are not sufficient to discriminate between all pairs of smooth 
billiard tables. Indeed, Lifshitz constructs (many) pairs of smooth domains (fR,^) which 
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have the same length spectra and the same wave invariants at corresponding pairs of closed 
billiard orbits 7 j of flj ( j = 1 , 2 ). We refer to [FKJ for proof that the domains with the same 
wave invariants in general have no eigenvalues in common. One can prove this by contra¬ 
diction using a deformation argument; see also EH] for the generalization to Schrodinger 
operators/ 

The domains are constructed in part from ellipses, which have very special billiard dy¬ 
namics, and are thus far from generic. They are also not real analytic. Thus they exemplify 
the strange behavior that can result in wildly non-generic cases mentioned in (D) of DO.21 It 
would be interesting to find spectral invariants that distinguish them. 

1.6. Residual spectral invariants and the Poisson-Wave trace. More generally one 
can consider zeta functions such as 

OO 

s{t, z) = Yl v v ' A '- ( Rez »°) ( 13 ) 

3 = 1 

and study the poles and residues of its analytic continuation in z E C. The residual spectral 
invariants are the residues at the poles. They have local geometric expressions, which relate 
spectrum and geometry. In particular the wave invariants above are residues of S(t,z ) at 
its poles in z for t 7 ^ 0 . Note that (J3]) is the boundary value of the holomorphic function 

S(t + ir) := e i{t+ir)Xg , r > 0 (14) 

j 

in the upper half plane {t + ir : r > 0}. More generally one may define the two-variable 
holomorphic function S(t + ir,z) = \J z e l ( t+lT ' >Xj , and study its (branched) meromorphic 
continuation from C + x C + (Rer > 0,Rez >> 0) to C x C. There may exist new spectral 
invariants from this analytic continuation. 

1.7. Non-residual or non-local spectral invariants. There exist many results on such 
non-local spectral invariants as the lowest eigenvalue Af or the Laplace determinant det A s . 
One can also consider special values of 0(f) or S(t,z ) at non-singular values. Non-local 
invariants are rarely computable in geometric terms. Examples which are the coefficients of 
positive integral powers of t in the heat trace Tre~ t/Xg at t — 0. I 11 this survey we concentrate 
on residual spectral invariants. 

2. Explicit formulae for wave invariants 

To explain what is meant by ‘local geometry’ we now give the formula for the ‘wave trace 
invariants’ or (equivalently) Birkhoff canonical form invariants around a closed geodesic 7 in 
the two-dimensional case. 

2.1. The principal wave trace invariant. A lot of information about (M,g) is already 
contained in the principal wave invariant fl 12 D . In |DG] (and its Appendix) it is proved 
that when Lsp(M , g) is simple, one may recover the lengths of all closed geodesics and the 
eigenvalues of all of the Poincare maps P 7 . This suggests already that the geodesic flows 
of isospectral Riemannian manifolds should not be far from symplectically equivalent in the 
sense discussed in 1 )0.21 at least in neighborhoods in S*M of closed geodesics. 
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2.2. Riemannian manifolds without boundary. The wave invariants at a closed geo¬ 
desic 7 are invariants of the germ of the metric at 7 . We introduce Fermi normal coordinates 
(s, y) along 7 and denote the corresponding metric coefficients by g t j . That is, the coordinates 
(s,y) correspond to exp 7 / s ) yv s where v s is the unit normal at the point 7 (s). 

We then consider orthogonal Jacobi fields along 7 , i.e. normal fields Y(s) = y(s)nu s where 
Y solves the Jacobi equation Y" + R(Y, Y)Y = 0 along 7 . The space of complex Jacobi fields 
along 7 is denoted The linear Poincare map is the monodromy map 

P 7 : Y(t) ->■ Y(t + L 7 ) (15) 

on this space. We refer to its eigenvectors Yj,Yj as Jacobi eigenvectors. 

We assume that the geodesic is non-degenerate in the sense that (det(J — P 7 ) 7 ^ 0). To 
simplify the exposition, we only consider the case where 7 is elliptic, i.e. where P 7 has unit 
modulus eigenvalues e iaj . We define the Floquet invariants 

Pi = (1 - e™*)- 1 . 

There is a similar story for hyperbolic closed geodesics. 

Theorem 2.1. [24] Let 7 be a strongly non-degenerate closed geodesic. Then the kth wave 
invariant at 7 has the form 

k / -Py;A:(®? d)ds 

' 7 

where: 

(i) I Ti k(s]g ) is a homogeneous polynomial of weight —k — 1 (under the scaling g —> e 2 g) in 
the data {y, y , D™ y g} with m = (mi,..., m n+ i) satisfying \m\ <2k + A; 

(ii) The degree of I r ^ in the Jacobi field Y = yv is at most 6 k + 6 / 

(in) At most 2k + 1 indefinite integrations over 7 occur in 1^.^; 

(iv) The degree of I Ti k in the Floquet invariants f3j is at most k + 2. 


The formula is simplest in dimension 2, where there is only one Floquet invariant (3. We 
use the notation r for the scalar curvature, t u for its unit normal derivative along 7 , t vv for 
the Hessian Hess(r)(z/, v). We denote by Y the unique normalized Jacobi eigenvector along 
7 and by Y its time-derivative. The subprincipal wave invariant a 7 0 is then given by: 

a 70 — 2'# [-®70;4(2/3 2 — (3 — —) + P 7 0 ; o] (16) 

where: 

(a) a 7i _i is the principal wave invariant (TT2]) : 

(b) L# is the primitive length of 7 ; a is its Morse index; P 7 is its Poincare map; 

(c) P 7 oy has the form: 


J pff 

B-yo-i = —~77 [ [a | Y\ 4 + b\ t\Y ■ Y | 2 + & 2 TRe (YY) 2 + c t 2 |Y | 4 + d t uv \Y\ 1 + e 5j 0 r\ds 

L* J Q 


— y 

L# ^ 


C v , 


0<m,n<3;m+n=3 


sin((n — m)a) 

|(1 — gi ( m -«)“)| 2 


»L# 


r I/ (s)Y m ■ Y n (s)ds 
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'l# 


E c * \mn 

0<ra,n<3;m+n=3 


Im 



r u (s)Y m • F n (s) 


r v (t)Y n ■ Y m (t)dt 


ds 


for various universal coefficients. The coefficients a, bi, b 2l c, d, Cj. mn are universal and can 
be determined from computable special cases. A somewhat computable case is that where 
7 is a rotationally invariant geodesic of length L of a surface of revolution, in which case 
t, 77 , t vu are constants and y = (we assume here that r = 1 on 7 ). 

If we iterate 7 (i.e. consider j k , traversing 7 k times), the Floquet invariants become 
independent and the two B l0 .j(j = 0,4) are independent. Thus each is a spectral invariant. 
This is a general phenomenon [G], Gf2l [Z3 ]. The coefficients B j0 .j are Birkhoff canonical form 
invariants. 


2.3. Bounded smooth plane domains. Euclidean plane domains hi are quite a different 
problem because the metric is flat (and known) and it is the boundary which is unknown. We 
express it locally as a graph y = f(x) over the x-axis. The role of closed geodesics is played 
by periodic trajectories of the biliiard flow on hi. In particular we consider ‘two-bounce’ 
or bouncing ball orbits, which are straight line segments hitting the boundary orthogonally 
at both endpoints. They always exist if the domain has an isometric involution or if it is 
convex. 

Rather than giving a formula for the subprincipal wave invariant, we give the top order 
terms (in derivatives of the boundary) for every wave invariant. Modulo terms involving 
< 2j — 2 derivatives, the wave trace invariants take the form 

B-yr,j-i = (4Lr )A r {ti)ii- l {2{wg2j,o) {h? r )if&\ 0) 

+4 K&‘i 3 ' o) W‘) j (/ (3) (°)/ 12j ~ 1) (°)) (17) 

{h£y-* E»Li('* 3 (/ (3) (0)/ (2j - i, (0))}. 

Here 

• h^, are the matrix elements of the inverse of the Hessian H 2r of the length function 
C in suitable (Cartesian graph) coordinates at 7 r ; 

• M r ( 0 ) is an ^-independent (non-zero) constant obtained from amplitude of the prin¬ 
cipal terms at the critical bouncing ball orbit. 

• w c 2 j,o (etc.) are certain non-zero combinatorial constants associated to certain Feyn- 

man diagrams (graphs). 

In this case, the Floquet exponents are buried in the Hessian matrix elements. A higher 
dimensional generalization of this formula is given in (HZ2j . 


3. Inverse results 

We now give an idea of how one can prove inverse spectral results using the wave trace 
invariants. 

3.1. Plane domains with an up-down symmetry. 
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3.1.1. Domains with one symmetry. The class consists of simply connected real-analytic 
plane domains hi satisfying: 

• (i) There exists an isometric involution a of Q; 

• (ii) a ‘reverses’ a non-degenerate bouncing ball orbit 7 -A y" 1 ; 

• (iii) The lengths 2 rL of all iterates y r have multiplicity one in Lsp(Q), and the 
eigenvalues of the linear Poincare map P 7 are not roots of unity; 

Let Spec(fi) denote the spectrum of the Laplacian of the domain with Dirichlet 
boundary conditions. 

Theorem 3.1. Spec: V l L ha M 1 ^ is 1-1. 

Here, Lsp{Vt) denotes the length spectrum, i.e. the set of lengths of closed trajectories 
of the billiard flow of H (see §2 for backgroun on billiards). By a bouncing ball orbit 7 is 
meant a 2-link periodic trajectory of the billiard flow. It corresponds to a line segment in 
the interior of D which intersects dfl orthogonally at both boundary points. By rotating and 
translating H we may assume that 7 is vertical and that A — (0, — ~). In a strip T e (AB), we 
may locally express dfl = r)0+ U dfl~ as the union of two graphs over the x-axis, namely 

dQ + = {y = f + (x), xe (-£,£)}, dQ~ = {y = f-(x), x e (-£,£)}. (18) 

The symmetry assumptions (i) - (ii) is that f+(x) = —/_(x). Condition (iii) on the multi¬ 
plicity of 2 L means that just the two orbits 7 , y -1 have length 2 L. In the elliptic case, its 
eigenvalues {e ±ia } are of modulus one and we are requiring that ^ ^ Q. I 11 the hyperbolic 
case, its eigenvalues {e ±A } are real and they are never roots of unity in the non-degenerate 
case. These are generic conditions in the class of analytic domains. 

Corollary 3.2. Let T> be the class of real analytic plane domains with an isometric invo¬ 
lution a satisfying: 

• (i) a reverses the shortest closed billiard trajectory 7 ; 

• (ii) rL 1 are of multiplicity one in LspifX). 

Then Spec: V ha M 1 ^ is 1-1. 

The proof of the corollary from the Theorem is just to observe that any shortest closed 
trajectory is automatically a bouncing ball orbit. Its length is a spectral invariant. 


3.1.2. Sketch of the proof. The key point is determine the 2j — 1st and 2 jth Taylor coefficients 
of the curvature at each reflection point from the j — 1 st wave trace invariant for 7 and its 
iterates . If the domain has the symmetries of an ellipse, one obtains the Taylor coefficients 
immediately from (TT71) . since the odd Taylor coefficients are zero. On the other hand, there 
is an obstruction to recovering the Taylor coefficients of / when there is only one symmetry: 
namely, we must recover two Taylor coefficients /^(O), /^ 2 jf- 1 ' ) (0) for each new value of j 
(the degree of the singularity). This is the principal obstacle to overcome. 

The expression inni for the Balian-Bloch invariants of 7 , 7 2 ,... consists of two types 
of terms, in terms of their dependence on the iterate r. They have a common factor of 
2rL(hHy~ 2 A r (0), and after factoring it out we obtain one term 


(y;) 2 {X 6f „)/ w, ( o) + 


( w a 2 j 7 1 ’ 3 - 0 ) 


2 . 1+1 


2 — 2 cosa /2 


/(3) (0)/( 2 ^—1) (0)} 
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which depends on the iterate r through the coefficient (hH) 2 , and one 

which depends on r through the cubic sums °f inverse Hessian matrix elements 

h 1 ^.- In order to ‘decouple’ the even and odd derivatives, it suffices to show that the functions 
(hH) 2 and EqLi(> 2 r ) 3 are, at least for ‘most’ Floquet angles a, linearly independent as 
functions of r E Z, i.e. that (h^) 2 is a non-constant function of r. It is 

convenient to use the parameter a = — 2 cos| and we write the dependence as /r 2 J r (a). 

We therefore define the ‘bad’ set of Floquet angles by 

2 r 

B = {a : the sequence {(hH(a))~ 2 ^(/r 2 )?(a)) 3 , ^ = 1,2,3,...} is constant in r}. (19) 

9=1 

Using facts about the finite Fourier transform and circulant matrices, we compute that 
B = {0,1, ±2}. For Floquet angles outside of £>, we can determine all Taylor coefficients 
f+\ 0 ) from the wave invariants and hence the analytic domain. 

3.2. A higher dimensional inverse result. It turns out that one can generalize the result 
to higher dimensions [HZ 2 j as long as the domain has Z 2 symmetries along all coordinate 
axes. The (Z/2Z) n symmetries are the maps 

CTj : (xi, ■ ■ ■, x n ) ->■ (xi,..., -Xj, Xj+i, • ■ ■, x n ). (20) 

The symmetry assumption implies that the intersections of the coordinate axes with 12 are 
projections of bouncing ball orbits preserved by the symmetries. 

We denote by Dl the class of all bounded real-analytic domains such that: 

' (i) cq : 12 —>■ 12 is an isometry for all j — 1 ,..., n; 

(ii) one of the coordinate axis bouncing ball orbits, called 7 , is of length 2 L 

(iii) the lengths 2 rL of all iterates 7 r (r = 1, 2, 3,...) have multiplicity one in Lsp(12); 

(iv) 7 is non-degenerate, i.e. 1 is not an eigenvalue of its Poincare map P 7 ; 
if 7 is elliptic and {e ±lQl , ...e ±4ari_1 } are the eigenvalues of P 7 , we 

further require that { 07 , ...,a n _i} are linearly independent over Q. We assume the 
_ same independence condition in the Hyperbolic case or mixed cases. 

( 21 ) 

Multiplicity one means that there exists precisely one closed billiard trajectory of the given 
length up to time reversal. Let Specs (12) denote the spectrum of the Laplacian of the 
domain 12 with boundary conditions B (Dirichlet or Neumann). 

Theorem 3.3. (HZ2j For Dirichlet (or Neumann) boundary conditions B, the map Specs '■ 
V L ->• is 1 - 1 . 

In other words, if two bounded real analytic domains 12i, 12 2 C M n possessing the symme¬ 
tries of an ellipsoid and satisfying the non-degeneracy and length assumptions of (| 2 Tjl have 
the same Dirichlet (resp. Neumann) spectra, then they are isometric. To our knowledge, the 
only prior positive result on the inverse spectral problem for higher dimensional bounded 
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domains is that a domain with the Dirichlet (or Neumann) spectrum of a ball must be a 
ball. 

The proof of Theorem 13.31 is similar in nature to the proof of Theorem 13.11 but uses the 
analysis in [ IHei| of the much more complicated dependence on Floquet eigenvalues in higher 
dimensions. 

3.3. Spectral determination of simple surfaces of revolution. A result which is some¬ 
what parallel to Theorem 13.31 is that convex analytic surfaces of revolution (satisfying a non¬ 
degeneracy hypothesis) are determined by their spectra (in this class). The profile curve of 
the surface of revolution is analogous to the ‘top half’ of the boundary of the plane domain, 
and the S 1 symmetry is analogous to the symmetry of the top and bottom halves. 

The precise class of metrics we consider are those metrics g on S 2 which belong to the 
class 1Z* of real analytic, rotationally invariant metrics on S 2 with simple length spectrum 
in the above sense and satisfying the following ‘simplicity’ condition 

• g = dr 2 + a(r) 2 d6 2 

• 3!ro : a'(r 0 ) = 0 ; 

• The Poincare map Vq of r = r$ is elliptic of twist type 

Convex analytic surfaces of revolution are examples, but there are others. The unique 
isolated closed geodesic (at distance r 0 ) is an elliptic orbit. 

Theorem 3.4. (|Z3] ) Suppose that g\ , g 2 are two real analytic metrics in the class 7 Z* on 
S 2 , i.e. ( S 2 , gi) are simple surfaces of revolution for which the meridian geodesic length is 
simple. Then Sp(A gi ) = Sp(A g2 ) implies gi = g 2 . 

We discuss the proof in 44.41 We note that the same kind of result should hold for 
Schrodinger operators (J2]) for fixed h — 1 , where one metric g € 1Z* is fixed and an S 1 - 
invariant potential V is varied. However, the ‘simplicity’ assumption probably needs to be 
placed on V and is probably similar to the monotonicity assumption in |DHV] , 

3.4. Isospectral deformations and spectral rigidity. As mentioned in the introduc¬ 
tion, one may at least heuristically consider the map Sp on the space of metrics, potentials, 
or domains. The derivative of Sp is the map 5g —> (<5Ai, 5X 2 , • • •) from the tangent space 
to metrics, potentials or domains to sequenes. One may consider variations of the spectral 
functions, wave invariants etc. Some early works on variations of eigenvalues are IS EOT. 

If 7 is an isolated, non-degenerate closed geodesic of g, then for any deformation g t of g, 
7 deforms smoothly as a closed geodesic of g t and one may define its variation 

^7 = fff \t=$L lt . ( 22 ) 

It is not hard to compute that 

L 1 = f gds, (23) 

J'y 

where g is viewed as a quadratic function on TM and 7 is viewed as the curve ( 7 (f), 7 ' (t)) 
in TM. 
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It follows that whenever the closed geodesics are non-degenerate and of multiplicity one 
in Lsp(M,g), an isospectral deformation preserves lengths and so 

gds = 0, V 7 . (24) 

In [GK] and in many subsequent articles, it is shown that there cannot exist such a deforma¬ 
tion among negatively curved surfaces. The proof is that there does not exist a (non-trivial) 
symmetric tensor g satisfying (|24l) . To our knowledge the analogous result for domains with 
boundary is unknown. The analogue of a negatively curved surface is a hyperbolic billiard 
table. Hyperbolicity refers to the Anosov property of the billiard flow. Such a table could be 
the complement of a disc (or more than one disc) in a negatively curved surface, or a plane 
domain whose boundary consists of a finite number of concave curves meeting at corners. 

Problem 2. Can one generalize the result of EH to hyperbolic billiards? Does a hyperbolic 
billiard table admit isospectral deformations? 



4. Geometric calculation of spectral invariants 

In this section, we briefly explain how one calculates the wave trace invariants explicitly 
in terms of geometric invariants or (in the boundary case) in terms of the boundary defining 
function. 


4.1. Parametrix ‘path integral’ versus canonical Hamiltonian approach. Roughly 

speaking, one has two approaches to constructing the ‘propagator’ U(t) = in quan¬ 

tum mechanics and studying its trace. One is the ‘Hamiltonian approach’ or ‘Birkhoff 
canonical forms’ approach in which one constructs the propagator locally around a closed 
geodesic 7 by adapting the canonical formalism (creation/annihilation operators, harmonic 
oscillators, etc.) to a tubular neighborhood of 7 . The canonical formalism is essentially the 
representation theory of the Heisenberg algebra and the symplectic algebra on L 2 (M n ). It is 
not obvious that one can ‘transplant’ the canonical formalism onto a Riemannian manifold, 
but this is what the theory of quantum Birkhoff normal forms achieves. Its origins lie in 
constructions of Gaussian beams by Babich, Lazutkin, Ralston and others (see |BB] ). Al¬ 
though Gaussian beams only exist for stable elliptic closed geodesics, the Birkhoff canonical 
form construction works locally around any non-degenerate closed geodesic (and even around 
degenerate ones). 

The second approach is the Lagrangian (or ‘path integral’) approach in which one con¬ 
structs a parametrix for the wave group. Unlike the canonical forms approach, there is noth¬ 
ing canonical about the parametrix construction. It is an art to construct a parametrix well 
adapted to a given problem. The parametrix approach is the traditional one and the reader 
is probably familiar with the heat kernel parametrix (if not with the term ‘parametrix’). 

Let us review the original wave kernel parametrix for (J7J) due to J. Hadamard (1921). On 
a Riemannian manifold without boundary, one can construct a Hadamard parametrix for 
cos,t\/~K{xpx) for small times, i.e. an approximation for the Schwartz kernel E(t,x,y) of the 
solution operator of the Cauchy problem, 


(f 2 -A)n = 0 
u\ t= o — f 


0 . ’ 


dt u \t=Q 


(25) 
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The parametrix is an oscillatory integral 


f 

Jo 


oo 


oo 


E(t,x,y ) ~ 


e' 


^j . 

,10(^2) m od G 00 


(26) 


i=o 


where the Hadamard-Riesz coefficients Uj are determined inductively by the transport equa¬ 
tions 



(27) 


The solutions are given by: 


U 0 (x,y) = Q 1/2 (x,y) 

U j+ i(x,y) = Q~ l/2 (x,y) f* s J 0(x, x s ) 1 / 2 A 2 U j (x, x s )ds 


(28) 


where x s is the geodesic from x to y parametrized proportionately to arc-length and where 
A 2 operates in the second variable. Also, r 2 (x, y) is the distance squared, dV g = 0(x, y)dy in 

"- 1 —j 

normal coordinates centered at x, and 0 r gg refers to the regularization of this distribution 
by analytic continuation (M. Riesz). 

The parametrix is a ‘geometric construction’ of the wave group and one may determine the 
singularities of the wave trace by expressing it in terms of the parametrix. The Hadamard 
parametrix above is useful to obtain the singularity at t = 0 and contains the information 
of the heat trace there. But the above parametrix is only valid for small t since it uses 
r 2 (x,y) and Q(x,y). To obtain information at much longer times such as t G Lsp(M,g ) 
one needs to construct a parametrix which is valid at that time around a closed geodesic. 
There are a number of clever ways to do this but none are simple except in special cases 
such as manifolds with no conjugate points. In that case one may follow the method of 
Selberg (1956) to construct the wave kernel on the universal cover and summing over the 
deck transformation group (see [Z1J for references and background). 

4.2. Domains with boundary. In the case of bounded domains, it is difficult to construct 
a parametrix that is valid in the glancing region of geodesics tangent to the boundary. A 
microlocal parametrix for the Dirichlet or Neumann even wave group E^(t) = cos ty/A b 
near a transversal reflecting ray was constructed by J. Chazarain EH (see also |GMl IPS] for 
more details ). 

To prove Theorem 13.11 the author used the construction of the wave group (or rather, 
Green’s function) using the method of layer potentials as in Balian-Bloch |BBll IBB2] rather 
than using a parametrix. This method is used universally in the physics literature to make 
calculations of eigenvalues and eigenfunctions and it is natural to use it for inverse spectral 
problems. 

We fix a non-degenerate bouncing ball orbit 7 of length L 7 , and let p G C^°(L y — e, L^ + e) 
be a cutoff, equal to one on an interval (L 7 — e/2, L 7 + e/2) which contains no other lengths 
in Lsp(G) occur in its support. The wave trace invariants associated to 7 are the coefficients 
of the asymptotic expansion as k —» 00 of 



( 29 ) 
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When 7,7 1 are the unique closed orbits of length L y , it follows from the Poisson relation 
for manifolds with boundary (see §3 and |GMj ) that (129|) has the expansion 

OO 

Trl n R p (k + hr) - e {ik ~ T)L ^ ^(S 7j + B y -ij)k- j , k -> oo (30) 

3 = 1 


with coefficients B y j, B y -ij determined by the jet of at the reflection points of 7 . The 
coefficients B y j, B y -i j are thus essentially the same as the wave trace coefficients at the 
singularity t — L y . 

As mentioned above, it is the art in inverse spectral theory to find a method by which 
the coefficients can be calculated explicitly. One might use a well-chosen parametrix around 
7 or use a canonical form. In [ Z4] , we use the layer potential method of [BBll lBB2j to 
express the Dirichlet (resp. Neumann) resolvent in terms of the ‘free resolvent’ Ro(k + ir) = 
— (A 0 + (k + zr) 2 ) -1 on R 2 , in the form, 

Rn(k + it) = Ro(k + ir) — V£(k + ir)(I + N(k + ir))~ lr yS£ tr (k + ir). (31) 

Here, 7 : H S (Q) —» H s ~ l ^ 2 {dVt) is the restriction to the boundary, and T>£(k + ir) (resp. 
S£{k + ir)) is the double (resp. single) layer potential, i.e. the operator from H s (dQ) —> 
H^ 1/2 (n) defined by 


Se(k + ir)f(x) = f m G 0 (k + iT,x,g)f(g)ds(g), 
V£(k + ir)f(x) = f m £^G 0 (k + ir, x, q)f(g)ds(g), 


(32) 


where ds(q) is the arc-length measure on dfl, where v is the interior unit normal to fl Also, 
S£ tr is its transpose (from the interior to the boundary), and Go(z,x,y) is the free Green’s 
function, i.e. the kernel of Ro(z). Further, 

N(k + ir)f(q) = 2 [ ^-G 0 (k + ir, q, q')f{q')ds(q') (33) 

Jon vvy 

is the boundary integral operator induced by VL One may simplify the trace by combining 
the interior and exterior problems. We write L 2 (R 2 ) = L 2 (Q) © L 2 (Q C ) and let R^, resp. 
(k+ir), R^ c (k+ir) denote the Neumann resolvent on the exterior domain, resp. the Dirichlet 
resolvent on the interior domain. We then regard Rp(k + ir) © R^(k + ir) as an operator 
on this space. We cycle around the layer potentials in (13TT) when taking the trace to obtain, 


Tr M 2 [R% c (k + ir) © R^(k + ir) - R 0 (k + ir)] = J- log det il + N(k + ir)J , (34) 

where the determinant is the usual Fredholm determinant. The main point then is to prove 
the existence of the following asymptotic expansion: 

Proposition 4.1. Suppose that L y is the only length in the support of p. Then, 


/ p{k — A) — log det(J + N{ A + ir))d\ ~ B y .jk J , 

Jm. dX . () 

where B r j are the wave invariants of 7. 
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One then shows that N(k + ir) is a semi-classical Fourier integral operator quantizing the 
billiard map on the boundary. Let us explain these terms. The billiard map (3 : B*dQ —> 
B*dQ on the ball bundle of dfl is defined as follows: given (y, rj) G B*dQ , with \r)\ < 1 we let 
(y, () G S*Q be the unique inward-pointing unit covector at y which projects to (y,rj) under 
the map Tg n Q —* T*dfl. Then we follow the geodesic (straight line) determined by (y, () to 
the first place it intersects the boundary again; let y' Gf denote this first intersection, and 
let rf be the projection of ( to B*,dQ. Then we define 

P(y,v) = 

We say that N(k + ir) quantizes the billiard map in that its Schwartz kernel is (almost) 
given by 

N(k + ir, y, y') ~ (k + ir)^ e^ k+iT ^ y - y %(k\y -y'\), 

where b{t) has an expansion in inverse powers of t as t —» oo, with leading term ~ t( _n+1 )/ 2 . 
This is proved in |HaZellZ4] , The kernel is a semi classical FIO of order 0 as long as \y — y'\ > 

_i 

k 2 log k. We separate out the tangential and transversal parts of N by introducing a cutoff 
of the form x(/c 1_<5 | q — q'\) to the diagonal, where 5 > 1/2 and where x *= Qj°(M) is a cutoff 
to a neighborhood of 0. We then put 


N(k + ir) = N 0 (k + ir) + Ni(k + ir), with 


(35) 


I N„(k + iT,q,c/) =x(k‘ ‘\q- ?'|) N(k + ir, q,q'), 

{ (36) 

N 1 (k+iT,q,q l ) = (1 - x(fc‘ ‘\q ~ q'\)) N(k + ir, q, q’). 

The term Aq is a semiclassical Fourier integral kernel quantizing the billiard map, while Nq 
behaves like an Airy operator close to the diagonal with the singularity of a homogeneous 
pseudodifferential operator on the diagonal. 

The advantage of this approach over the canonical forms method or using a Chazarain (or 
other parametrix) is that the kernel N(k, q, q') is independent of the boundary except for the 
restriction of the layer potential to the given boundary. In some sense the expansion above 
is ‘canonical’. The phase \y — y'\ on dkl x dkl is the generating function of the billiard map as 
long as is convex in the sense of symplectic geometry. If fails to be convex, it generates 
a canonical relation which is larger than the graph of (3, causing some technical problems 
that can be handled [HaZej . One also must take into account the part of the kernel very 
close to the diagonal, which is particularly important if one is studying billiard trajectories 
that creep along the boundary but it is not important if one considers trajectories that are 
transverse to the boundary such as bouncing ball orbits. 

At least formally, we then expand (/ + N(k + ir)) -1 in a finite geometric series plus 
remainder: 


M 0 

(I+N(X+ir)) -1 = J](-1) M N(X) M + (-1) M ° +1 N(X) Mo+1 (I+N(X+ir)) -1 . (37) 

m =o 


(N 0 + N,) M 


° 

cr :{ l ,..., M }—>{ 0 , 1 } 


o • • • O N a{ M) ■ 


We then write 


(38) 
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We regularize N M by eliminating the factors of Nq from each of these terms. This is obviously 
not possible for the term N ( f f but it is possible for the other terms. By explicitly writing 
out the composition in terms of Hankel functions and using the basic identities for these 
special functions, we prove that N 0 oNi o xo(k + ir, tpi, tp 2 )) is a semiclassical Fourier integral 
operator on dfl of order —1 associated to the billiard map. Thus, composition with N 0 lowers 
the order. The remaining terms N ^ 1 , when composed with a cut-off to 7 , do not contribute 
asymptotically to the trace. 

The successive removal of the factors of N 0 thus gives a semi-classical quantization of 
the billiard map near 7 . We then calculate the traces of each term by the stationary phase 
method to obtain the expansion (15UD . The coefficients containing the top derivatives for each 
power of k are rather straightforward to calculate explicitly. 


4.3. Birkhoff canonical form. So far we have described parametrix methods for evalu¬ 
ating wave trace coefficients. We now describe the Birkhoff canonical form and the sense in 
which it is canonical. Let 7 be a non-degenerate closed geodesic on an n -dimensional Rie- 
mannian manifold, and at first let us assume it to be elliptic. For each transversal quantum 
number q G Z n_1 , there was an approximate eigenvalue of the form 


^kq — Tkq "P 


Pl(g) 

Tkq 


+ 


P 2 (q) 


+ ... 


kq 


(39) 


where 

1 n 1 

r kq = ^( 27 Tk + + 2 )^) 

3 =1 

where the coefficients are polynomials of specified degrees and parities. We refer to [BB ] for 
a clear exposition and for details. 

A natural question is whether the wave invariants a 7J - of (fTOj) can be determined from the 
quasi-eigenvalues (l39]h i.e. from the polynomials Pj(q). The answer is ‘yes’ and in effect (1391) 
is the canonical form. 

To put A into normal form, is first to conjugate it into a distinguished maximal abelian 
algebra A of pseudodifferential operators on a model space, the cylinder S\ x M”, where S\ 
is the circle of length L. The algebra is generated by the tangential operator D s := 4^ on 
together with the transverse action operators. The nature of these action operators depends 
on 7 . When 7 is elliptic, the action operators are harmonic oscillators 

Ii = Ii(y,D v ) 


while in the real hyperbolic case they have the form 


I3 — VjDyj + Dy^t/j. 

When 7 is non-degenerate, they involve some mixture of these operators (and also complex 
hyperbolic actions) according to the spectral decomposition of P 7 . For notational simplicity 
we restrict to the elliptic case and put 

1 n 

H a ■ 7 ^ ^ ^ 

z k =1 

where e ±lctk are the eigenvalues of the Poincare map P 7 . 
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To put A into normal form is to conjugate it to the model space and algebra as a function 
of D s and the action operators. The conjugation is only defined in a neighborhood of 7 in 
T*M — 0, i.e. one constructs a microlocally elliptic Fourier Integral operator W from the 
conic neighborhood of M + 7 in T*N 7 — 0 to a conic neighborhood of TfSf in T*(Sf x R n ) 
such that: 


Ws/KpW~ l = v 


Pi{h 


LV 


+ 


P2(Il} ■ ■ ■ ) In') Pk+l(Ili ■ ■ ■ ) In) 


(Lvy 


(LV) 


fc+i 


(40) 


where the numerators Pj(Ii ,..., I n ),pj(ii ,..., I n ) are polynomials of degree j+1 in the variables 
Ii,...,I n , where W^ 1 denotes a microlocal inverse to W. Here, V = D s + \H a . The kth 
remainder term lies in the space (B*to 02 (k+ 2 -j)'^ 1 ~P where denotes the pseudo-differential 
operators on the model space of order s and where C> 2 (k+ 2 -j) denotes the operators whose 
symbols vanish to the order 2(k + 2 — j) along 7 . We observe that (j40l) is an operator version 
of (BSD. 


The inverse result of [G], see also [Z3l [Z4]] is: 


Theorem Let 7 be a non-degenerate closed geodesic. Then the quantum Birkhoff normal 
form around 7 is a spectral invariant; in particular the classical Birhoff normal form is a 
spectral invariant. 

In other words, one can determine the polynomials p 3 (I \, ...,/„) from the wave trace in¬ 
variants of A at 7 . The key point is that one can construct the normal form explicity by 
re-scaling the metric around 7 and making a step-by-step conjugating to Hermite operators 
and creation/annihilation operators first define by Babich, Lazutkin and others (see [BB] 
for background). Thus, the algorithm is quite concrete and leads to computable spectra 
invariants of Theorem 12.11 But to date, the formulae have not been applied to any concrete 
inverse spectral problems due to their complexity. In [Z21IHZ2] we found that inverse results 
could be proved in the boundary case by just keeping track of the highest order derivative 
terms in the wave invariants. It is possible that the argument could be generalized to some 
metric cases. 

4.4. Proof of Theorem 13.41 on spectral determination of analytic simple surfaces 
of revolution. As mentioned above, surfaces of revolution are the surfaces most analogous 
to bounded plane domains in that both are defined by a profile function of one real variable. 
This suggests that one should be able to prove an analogue of Theorem 13.11 for analytic 
surfaces of revolution. In fact, it was proved first and by a normal forms method. To our 
knowledge, it remains the only inverse spectral result for metrics which was proved by using 
Birkhoff normal forms, and therefore we briefly recall the method. Recently, similar methods 
have been used in the K potential problem for (J2J) in several settings [CdVGl ICdV 2 [ IGW . GUj. 

The proof is based on quantum Birkhoff normal forms for the Laplacian A. But the special 
feature of simple surfaces of revolution is that there exists a global Birkhoff normal form as 
well as local ones around the critical closed orbit or the invariant tori. It was constructed 
in |CdV4] and expresses the fact tht A is a toric integrable Laplacian in the following sense: 
there exist commuting first order pseudo-differential operators A, G such that: 
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• the joint spectrum is integral, i.e. Sp(£) C Z 2 flT + {/i} where T is the cone J 2 > |/i| 
in M 2 . 

• The square root of A is a first order polyhomogenous function \fK = H{I \, If) of the 
action operators. 

By polyhomogeneous, we mean that H has an asymptotic expansion in homogeneous 
functions of the form: 

H + + Hj (■ rl ) = r 3 Hj (/). 

The principal symbols Ij of the Ifs generate a classical Hamiltonian torus action on T*S 2 — 0. 
There is no term of order zero. It follows that 

Sp(y/K g ) = {H{N + n ): N e z 2 n r 0 }, 

where the eigenvalues have expansions 

A n ~ Hi(N + p) + H_i(N + //) + .... 

Theorem 4.2. f|Z2jj Tet (S 2 ,g) be an analytic simple surface of revolution with simple 
length spectrum. Then the normal form 2 ) is a spectral invariant. 


The normal form and the proof are very different from the non-degenerate case in K31 
E3 ESJ. Recently, similar kinds of results have been proved in other integrable settings 

Icpyifcwj 

To complete the proof of Theorem 13.4[ we need to show that H determines a metric in 
1Z. As in the bounded domain case outline above, we need to calculate the normal form 
invariants. It turns out to be sufficient to calculate Hi = H and iT_ 1 in terms of the metric 
(i.e. in terms of a(r)) and then to invert the expressions to determine a(r). In [Z2] H and H_ 1 
were calculaed using by studying asymptotics \fK = H(h, I 2 ) along ‘rays of representations’ 
of the quantum torus action, i.e. along multiples of a given lattice point (■ n 0 ,k 0 ). In effect, 
the no parameter is like hr 1 and we reduce to a one-dimensional Sturm-Liouvillc problem. 
Thus, Theorem 14.21 converts the hard problem into the easier problem of determining the 
potential for semi-classical Schrodinger operators. In a rather standard way, one finds that 

f (E — x)\dp,(x) 

Jr 


is a spectral invariant, where /j is the distribution function p,{x) := |{r : < x}| of 

with | • | the Lebesgue measure. This Abel transform is invertible and hence 


dp(x) 



d 1 
dr a(r) 2 


l dx 


and therefore 


J(x) := 


E 

r:a(r)=x 


1 

a'(r )| 


are spectral invariants. By the simplicity assumption on a, there are just two solutions of 
a{r ) = X] the smaller will be written r_(x) and the larger, r + (x). Thus, the function 


J(x) 


1 

a'(r-(x)) | 


+ 


1 

a'(r + (x))| 


(6.9) 
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is a spectral invariant. 

By studying if_i, we find in a somewhat similar way that 

K{x) = |a , (r_(x))| + |a , (r + (x))| 

is a spectral invariant. It follows that we can determine a'(r + ±(x)). 

Since both metrics g\ and r/ 2 are assumed to belong to 1Z*, they are determined by their 
respective functions a 3 {r). We conclude that a\ = a 2 and hence that (]\ = <? 2 . 

In more recent work on the inverse problem for the h-spectra of Schrodinger operators, 
the methods are similar to the final step in the proof above but allow potentials which have 
more critical points. It is then a challenge to recover the potential from its distribution 
function and the higher invariants. Vice-versa, the methods might be used for more general 
inverse problems for metrics on surfaces of revolution in which one relaxes the simplicity 
assumption. 


5. ISOSPECTRAL DEFORMATIONS OF DOMAINS 


An isospectral deformation of a plane domain Q 0 is a one-parameter family hL of plane 
domains for which the spectrum of the Euclidean Dirichlet (or Neumann) Laplacian A e is 
constant (including multiplicities). We say that is a C 1 curve of C°° plane domains if 
there exists a C 1 curve of diffeomorphisms ip £ of a neighborhood of O 0 C M 2 with p 0 = id and 
with = 93 £ (h2 0 ). The infinitesimal generator X = is a vector field in a neighborhood 
of Q 0 which restricts to a vector field along <9f2 0 ; we denote by X v = pv its outer normal 
component. With no essential loss of generality we may assume that <^ e |an 0 is a map of the 
form 

x G <9f2o —y x + p £ (x)is x , (41) 

where p £ G C^QO, £ 0 ], C°°(dQ 0 )), > 0 and p 0 = 0. We put p(x) = Sp(x) := j- £ \ £=Q p £ {x). 

An isospectral deformation is said to be trivial if Vt £ ~ Oo (up to isometry) for sufficiently 
small e. A domain Oo is said to be spectrally rigid if all isospectral deformations are 
trivial. The domain O 0 is called infinitesimally spectrally rigid if p — 0 (up to rigid motions) 
for all isospectral deformations. 

The variations of the eigenvalues are given by Hadamard’s variational formulae, 

{ IdQ(o) P \du^j(0)\ d fi\ 2 dA, Dirichlet 

(42) 

foa 9 | </b (0) l£>f 2 (o) | 2 dA, Neumann. 

Hence, the infinitesimal deformation condition is that the right hand sides are zero for all 
j. To normalize the problem, we assume with no loss of generality that the deformation is 
volume preserving, which implies that 


pdA = 0. 


(43) 


'on 


Any such p defines a volume preserving deformation of f2. 

Thus, the infinitesimal deformation is orthogonal to all boundary traces of eigenfunctions: 


Ai(0) = OVj, 


f d Q ( o) P \d u pj(0)\an\ 2 dA = 0, Dirichlet 
Sxi P l^i(0)|an(o)| 2 dA = 0, Neumann. 


(44) 
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We may rewrite these conditions in terms of the boundary values of the wave kernel: 


as saying that 


E b (t , x, x ) 


d Ux d Uy U(t,x,x)\ XIEdn , Dirichlet 
U(t,x,x)\ X £dn, Neumann 



E b (t , x, x)p(x)dA(x) 


0 ,Vt. 


( 45 ) 


(46) 


5.1. Spectral rigidity of the ellipse. 

2 2 

Problem 3. /s an ellipse E a j, given by ^ = 1 determined by its Dirichlet or Neumann 

spectrum? Is it spectrally rigid, i.e. does there exist a non-trivial isospectral deformation fl t 
of E a ^ ? Since the area and perimeter are spectral invariants it is understood that fl t has 
fixed area and perimeter. 


In the case of circles, the answer was shown to be ‘yes’ by M. Kac, since the disc is 
extremal for the isoperimetric inequality. But the answer remains unknown for general 
ellipses. Ellipses are special because they have completely integrable billiard flows, and 
(since G. D. Birkhoff) are conjectured to be the only smooth plane domains with integrable 

billiards. 

In [HZ], H. Hezari and the author proved a partial rigidity result which breaks out of the 
class of real analytic domains, but does assume that the domains all possess the symmetries 
of the ellipse. Namely, ellipses are spectrally rigid among C°° plane domains with the 
symmetries of an ellipse. Here ‘rigidity’ means that there do not exist smooth curves fit 
in the family of isometry classes of smooth domains with the spectrum of the ellipse, at 
least if the curves are ‘non-flat’ in that the Taylor expansion of the domain at the endpoins 
of the axes is non-zero in the deformation variable t. This non-flatness assumption is an 
unexpected and annoying technicality which it would be nice to remove. The result is 

Theorem 5.1. Suppose that f2 0 is an ellipse, and that e —>■ is a C°° Dirichlet (or 

Neumann) isospectral deformation through Z 2 x Z 2 symmetric C°° domains. Then p e must 
be flat at e = 0. In particular, there exist no non-trivial real analytic curves e —>■ of 

Z 2 x Z 2 symmetric C°° domains with the spectrum of an ellipse. 

We used the Hadamard variational formula in the proof. By taking the variation of the 
wave trace, we proved: 

Theorem 5.2. Let Qo C W 1 be a C°° convex Euclidean domain with the property that the 
fixed point sets of the billiard map are clean. Then, for any C 1 variation of through C°° 
domains fl E , the variation of the wave traces 5Tr cos (t\/—A), with Dirichlet (or Neumann) 
boundary conditions is a classical co-normal distribution for t m|<9f2 0 | (m & Z) with 
singularities contained in Lsp(flo)- For each T e Lsp(flo) for which the set Ft of periodic 
points of the billiard map j3 of length T is a d-dimensional clean fixed point set consisting 
of transverse reflecting rays, there exist non-zero constants Cp independent of p such that, 
near T, the leading order singularity is 

5 Tr cos (ty/~ A) ~ ^ Re { ( ^ C r f p 71 dp r ) (t — T + f0 + )' 2- ^}, 

rcF T 






THE INVERSE SPECTRAL PROBLEM 


23 


modulo lower order singularities. The sum is over the connected components T of Ft- Here 

6 - iU=o and 7i(9, C) = \/l - Kl 2 - 

Here 71 is a certain computable function on B*dfl. 

the billiard flow and billiard map of the ellipse are completely integrable. In particular, 
except for certain exceptional trajectories, the periodic points of period T form a Lagrangian 
tori in S*VLq. The exceptions are the two bouncing ball orbits through the major/minor axes 
and the trajectories which intersect the foci or glide along the boundary. The fixed point 
sets of <1> T intersect the co-ball bundle B*dQ 0 of the boundary in the fixed point sets of the 
billiard map (3 : B*dfl 0 —» B*dfl 0 (for background we refer to jPSl IGM1 HZ for instance). 
Except for the exceptional orbits, the fixed point sets are real analytic curves. For the 
bouncing ball rays, the associated fixed point sets are non-degenerate fixed points of (3. 

The densities dp r on the fixed point sets of (3 and its powers are very similar to the 
canonical densities defined in Lemma 4.2 of |DG] . and further discussed in jGMl PT], The 
constants Cp are explicit and depend on the boundary conditions. Also, Ft is the curve 
of periodic points of the billiard map /3 on B*dQ of length T, where as above the length 
is defined by the length of the corresponding billiard trajectory. The dimension d is the 
dimension of F T . In the case of the ellipse, for instance, d — 1; the periodic points of a given 
length form invariant curves for (3. 

Theorem 5.3. Suppose that is an ellipse, and that Q, e is a C 1 Dirichlet (or Neumann) 
isospectral deformation of f2 0 through C°° domains with Z 2 x Z 2 symmetry. I. e. is invariant 
under (x , y) —>■ (±x, ±y). Thenp = 0. 

This implies that ellipses admit no isospectral deformations for which the Taylor expansion 

_ i 

of p e at £ = 0 is non-trivial. A function such as e ? for which the Taylor series at e = 0 
vanishes is called ‘flat’ at e — 0. 

Since the final step of the proof uses results of ( GM] . we briefly review the description of 

2 2 

the billiard map of the ellipse flo := 7 - + \ — 1 (with a > b > 0) in that article. I 11 the 
interior, there exist for each 0 < Z < b a caustic set given by a confocal ellipse 

x 2 y 2 

-h — = 1 

E + Z Z 

where E = a — b, or for — E < Z < 0 by a confocal hyperbola. Let (g, f) be in B*dflo and let 
(g, £) in be the unique inward unit normal to boundary that projects to (g, (). The line 
segment (g, r() will be tangent to a unique confocal ellipse or hyperbola (unless it intersects 
the foci). We then define the function Z(q,f) on B*dVt 0 to be the corresponding Z. Then 
Z is a /3-invariant function and its level sets {Z = c} are the invariant curves of /?. The 
invariant Leray form on the level set is denoted duz (see [GM] . (2.17), i.e. the symplectic 
form of B*dflo is dq A df = dZ A duz- A level set has a rotation number and the periodic 
points live in the level sets with rational rotation number. As it is explained in [ GM] (page 
143) the Leray form duz restricted to a connected component T of Ft is a constant multiple 
of the canonical density dp r- 

As mentioned in the introduction, the well-known obstruction to using trace formula 
calculations such as in Proposition 15.21 is multiplicity in the length spectrum, i.e. existence 
of several connected components of Ft- A higher dimensional component is not itself a 
problem, but there could exist cancellations among terms coming from components with 
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different Morse indices, since the coefficients Cp are complex. This problem arose earlier in 
the spectral theory of the ellipse in [GMj . Their key Proposition 4.3 shows that there is a 
sufficiently large set of lengths T for which F T has one component np to (q,C) (<Z, —C) 

symmetry. Since it is crucial here as well, we state the relevant part: 

Proposition 5.4. (see [GMj . Proposition f.3): Let Tq = |<9f2o|- Then for every interval 
(roT 0 — e, mT Q ) for m = 1,2,3,... there exist infinitely many periods T e Lsp(G 0 ) for which 
F t is the union of two invariant curves which are mapped to each other by (q, () —>■ ( q , —(). 

Since for an isospectral deformation 5 Tr cos (t\/—A) = 0, we obtain from Proposition 15.21 
the following 

Corollary 5.5. Suppose we have an isospectral deformation of an ellipse G 0 with velocity 
p. Then for each T in Proposition \5.f\ for which Ft is the union of two invariant curves Ti 
and r 2 which are mapped to each other by (g,C) —>■ (g, —C) we have 

/ P 7i du z = 0, J = 1,2. 

Jrj 

Proof. From Proposition 15.21 we get 

2 f 

Re { ( C Tj / p 71 dpr :i ) (t - T + i0)~ 2 ~ a } = 0 . 
j =1 

Since p and 71 are invariant under the time reversal map (q,C) ~t (q, —(), the two integrals 
are identical. Also by directly looking at the stationary phase calculations it can be shown 
that the Maslov coefficients Cpi and Cp 2 are also the same. Thus the corollary follows. □ 


5.2. Abel transform. The remainder of the proof of Theorem 15.31 is identical to that of 
Theorem 4.5 of |GM] (see also Ell). For the sake of completeness, we sketch the proof. 

Proposition 5.6. The only Z 2 xZ 2 invariant function p satisfying the equations of Corollary 
\5J\isp = 0 . 


Proof. First, we may assume p — 0 at the endpoints of the major/minor axes, since the 
deformation preserves the Z 2 x Z 2 symmetry and we may assume that the deformed bouncing 
ball orbits will not move and are aligned with the original ones. Thus p(±\/a) = p(±Vb) = 0 . 

The Leray measure may be explicitly evaluated (see 2.18 in [ GMj ). By a change of variables 
with Jacobian J, and using the symmetric properties of p, the integrals become 


A(Z) = 



p{t) 7 ! J{t)dt 
yt-(b-z) 


(47) 


for an infinite sequence of Z accumulating at b. The function A(Z) is smooth in Z for Z 
near b. It vanishes infinitely often in each interval (6 — e, 6 ), hence is flat at b. The kth Taylor 
coefficient at b is a 

A^ k \b) = J p(t) 7l J{t)r k ~\dt = 0 . (48) 

Since the functions t~ k span a dense subset of C[b, a], it follows that p = 0 . 


□ 


We refer to 1 17.21 for some further problems regarding the isospectral determinantion prob¬ 
lem for ellipses. 
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5.2.1. Spectral rigidity of convex analytic domain with the symmetries of an ellipse? We 
observe that the proof of Theorem 15.11 only uses the hypothesis that the domain is an ellipse 
in two steps. First, in Theorem 15.21 it uses the very special kind of singularities of the 
wave trace due to the foliation of B*dE by invariant curves for the billiard flow. Then in 
Proposition 15.61 it uses the explicit formula (14711 for the integrals over the invariant curves. 
In the case of general convex smooth domains, there exist invariant sets for the billiard map 
due to periodic creeping orbits along the boundary, but they are much more complicated 
than in the completely integrable case of the ellipse. Still, it may be possible to analyze the 
analogous integrals A(Z). 

6. Inverse scattering, phase shifts, resonance poles 

In this section we consider the inverse phase shift problem for semi-classical Schrodinger 
operators ([2]) on (with the Euclidean Laplacian). The scattering matrix 

S h : L 2 (§ d_1 ) -A L 2 ^- 1 ) 

can be defined in terms of generalized eigenfunctions as follows. For ipi n E C ,0O (§ d ^ 1 ), there 
is a unique solution to Hu = 0 satisfying 

u = (e-^r/V in (w) + e ^/V 0Ut (_u;)) + 0 (r~^)A), ( 49 ) 

By definition 

S h {pin) ■= d <d - x)l \ ou t■ (50) 

Below, we will refer to pi n and ip ou t as the incoming and outgoing data of u. It is not 
hard to show that 

S h y{E) = 5^(1), (51) 

where h = h/\/E, V = V/E, and Sh'y'(E') denotes the scattering matrix for (h') 2 A+V — E’. 
Using (ISTjh we may assume that E — 1 . 

Phase shifts are the eigenvalues for the eigenvalue problem 

S h <Ph,n = e i0h ' n Vh,n- (52) 

Problem 4. Can one determine the potential V (in some reasonable class) from the phase 
shifts {e l/3h ’ n }(f =1 for fixed energy but as h varies? 

The simplest problem is to assume maximal symmetry, i.e. that the potential is radial. 
This effectively reduces the problem to one dimension. The scattering matrix then commutes 
with the action of SO(d) and A S d-i and thus has just one eigenvalue in each eigenspace of 
the Laplacian. When d = 3, the eigenvalues of A 52 on the sphere are known as the total 
angular momenta, and thus the phase shifts are traditionally denoted 5^(k) where the angular 
momentum is £{£ + 1 ) and k = hr 1 . 

The question of finding a radial potential V from its phase shifts has been much stud¬ 
ied since the 1960's. The phase shifts are closely related to the Jost functions fe(k) = 
\fe(k)\e~ lSe ^\ The first issue is whether ([2]) has any L 2 -eigenfunctions or ‘bound states’. In 
classical scattering the corresponding issue is whether there exist trapped Hamilton trajecto¬ 
ries or whether they all escape to infinity. It one knows all of the eigenvalues E n ( K) and the 
phase shifts 5i(k) for all k > 0 then one can construct \fg(k)\. Equivalently one knows the 
ratio fe(—k)/fe(k). The problem of determining fe(k) from the ratios is a Riemann-Hilbert 
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problem. Using the Gelfand-Levitan solution of the inverse spectral problem in dimension 
one, it has been shown that if there are N eigenvalues then there exists an N-parameter 
family of potentials with the same phase shifts as k varies and £ is fixed and same eigen¬ 
values [Nj IN21 |N3j . In particular, the potential is uniquely determined if (J2J) has no L 2 
eigenfunctions. If instead one fixes k and lets £ vary, one can also solve the problem but not 
uniquely. 

To our knowledge, the problem of determining V from the /3h, n in non-radial cases is open 
and the analogues of the methods described above for Laplacians on compact manifolds have 
not been discussed. 

A closely related problem is to study the Laplacian A of the Euclidean metric in the 
exterior of a convex obstacle O. Thus, we let 0 = and consider the scattering matrix 

Sh for h 2 A on L 2 (U) with Dirichlet or Neumann boundary conditions. 


Problem 5. Can one determine 0 from the fixed energy phase shifts e l ^ h ’ n ?. 


The inverse problem for obstacles is trivial in the radial case but very non-trivial in general. 
There exist results on determining the potential from the scattering kernel (Theorem 2 of 
|Majl9 76, 1G3] ). but to our knowledge the phase shift problem has not been studied. The 
method of [G3] is to determine the so-called sojourn times from the full scattering kernel 
Sh(aj,9) as h —* 0 and to show that it dtermines the Gaussian curvature of the surface. In 
the case of a convex obstacle, the curvature function determines the obstacle, i.e. the shape 
of the scatterer. The question of how much of the shape of the scatterer one might recover 
from the phase shifts alone does not seem to have been studied in the non-radial case. 

As J. Gell-Redman emphasized to the author, the so-called interior-exterior duality estab¬ 
lishes a relation between the interior Dirichlet eigenvalues and the exterior Dirichlet phase 
shifts [ EPlJ . From the phase shifts one can determine the interior Dirichlet eigenvalues. Hence 
positive results on the inverse problem for domains imply positive results for phase shifts. 

We refer to |Maj 1976] , IChSl ICT1IKKS1IMT1 IM IN21IN31 ISabj for background on inverse re¬ 
sults in dimensions > 1 and to ( GHZ . PG11112013 for some recent results on the distribution 
of phase shifts. Some further references in dimension one are given in the introduction. 


7. Problems and partial solutions 

In this section we state some problems which are at least partly open. They are intended 
both to give a sense of the current state of the art and also to suggest problems on which it 
seems likely one can make progress. 

The following is probably the most natural: 


Problem 6. Is a convex analytic domain spectrally rigid in the class CA of convex smooth 
domains? More difficult: Is a C°° convex domain spectrally rigid? 


There are several possible approaches. One is to develop the proof of Theorem 15.11 Only 
in the last step is the geometry of the ellipse used. It will change radically for a general 
convex domain. In particular one must use the invariant sets which arise in general. 
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In view of the number of inverse spectral results which use the assumption that the domains 
are real analytic, we pose: 

Problem 7. Suppose H] and fi 2 ore isospectral plane domains. If 01 is real analytic, is fi 2 
real analytic? Is this true if we assume both domains have the symmetry of an ellipse? 


The reason for posing the problem is that, when a domain has the symmetry of an ellipse, 
one can determine the Taylor expansion of the domain at the endpoints of the two axes 
from the spectrum. Hence if the expansion converges for the first domain, the Taylor series 
converges for the second domain. But it might not actually equal the second domain. More 
precisely, we represent the top and bottom of the domains as graphs y = f±(x) over one of 
the axes and then the Taylor expansion of f±(x) j — 1,2 is determined by the spectrum. 
Under the symmetry assumption f 3 _ = — f 3 + . For instance, in the case of an ellipse E a b, 

f+(x ) = byjl — We note that the radius of convergence of f± at x — 0 is |x| < a, i.e. 
until the graph ceases to be a graph by turning vertical. The length of the axes is a spectral 
invariant and therefore it must coincide with the first domain. 

The question arises whether fi 2 is locally defined by f±. It is possible that f± — f± + if± 
where if vanishes to infinite order at the endpoints of the bouncing ball orbit. The wave 
invariants still coincide but it is not clear that the wave trace itself does. One might need 
“exponentially small corrections” to the wave trace to determine that, or else a proof that the 
wave trace expansion converges. It is not clear that this is possible for isospectral surfaces. 
A simpler problem is 

Problem 8 . Suppose Hi and fi 2 are isospectral plane domains. If Hi is real analytic, is fi 2 
locally real analytic at the endpoints of a bouncing ball orbit? Is this true if we assume both 
domains have the symmetry of an ellipse? 


If so, we then ask: Suppose that fi 2 is isospectral to a real analytic Hi and that <9H 2 = <9fii 
on some interval x E (—e,e) around the endpoints of a bouncing ball orbit. Must then 
fi 2 = Hi? 

Another question concerns the normal form in the analytic case. 


Problem 9. Does the quantum Birkhoff normal form at a hyperbolic geodesic 7 converge if 
the metric is real analytic? Does the conjugating map (unitary operator) converge? 


It is proved in [Mo] that this is true for the classical Birkhoff normal forms in dimension 2. 
Some hints that the answer is positive for quantum normal forms in dimension 2 are given 
in Sj, [Aj ISa.. PM] IT], If the answer is ‘yes’ then the manifolds are locally Fourier-isospectral 
near 7 . In particular the geodesic flows are locally symplectically equivalent. 


Problem 10. Do isospectral plane domains have symplectically equivalent billiard maps? If 
a domain H has the same eigenvalues as the ellipse E a b , does it have a completely integrable 
billiard map? 


Results of Siburg jS] and of the author and G. Forni show that the billiard flow would be 
C° conjugate to that of the ellipse, at least if one has an isospectral deformation. But there 
is a long distance between C° and smooth conjugacy. 








28 


STEVE ZELDITCH 


7.1. Problems on surfaces. 

Problem 11. If g is a metric on S 2 such that the distinct eigenavalues X k of A g have the 
same multiplicities m(X* k ) = 2k + 1 as the standard metric g 0 , is g = g 0 ? 

7.2. Questions about oscillatory integrals. 

Problem 12. If an oscillatory integral f Rn ae lktf dx whose phase has a single critical point at 
x = 0 has the stationary phase asymptotics of a Morse function with a single critical point, 
is it a Morse function? Similarly, if it has the stationary phase asymptotics of a Bott-Morse 
phase function with a non-degenerate critical manifold, is it Bott-Morse? 

The reason for posing this problem is that wave trace invariants are obtained by applying 
stationary phase to a parametrix for the wave group, and one only knows the wave trace 
invariants explicitly if the closed geodesics are non-degenerate (i.e. if the phase of the integral 
is Morse) or ‘clean’ (if the phase is Bott-Morse). If a second domain or metric has the same 
wave trace invariants as the first, then its wave trace expansion at a closed geodesic is the 
same as one for first, given by stationary phase for a Morse function. It would simplify the 
inverse problem considerably if we could know that that the closed geodesic of the unknown 
metric or domain was also Morse (or Bott-Morse). For instance in the case of the ellipse, 
the phase functions are Bott-Morse. 

As a potential application, we pose: 

Problem 13. Suppose that hi is isospectral to an ellipse. Do the periodic orbits of the 
billiard map j3 on B*dLl form smooth invariant curves? Similarly, if ( S 2 ,g ) is a metric 
on S 2 which is isospectral to a convex surface of revolution, do the periodic geodesics form 
smooth invariant Lagrangian tori? 

The point is that the wave trace invariants are the same as for the ellipse, resp. ( S 2 ,g ) 
(here we ignore issues of length spectral multiplicity). As we see from Theorem 15.21 the wave 
trace has the singularities of a Bott-Morse oscillatory integral with clean fixed point sets of 
dimension 2 (in S*M ) resp. 1 (in B*dD). A positive answer to the Problem about oscillatory 
integrals with the stationary phase asymptotics of a Bott-Morse function would imply that 
the phase has to be Bott-Morse for the unknown metric or boundary and therefore would 
be smooth invariant submanifolds. 

A related more ambitious question is: 

Problem 14. Suppose that is isospectral to an ellipse. Is B*dLl foliated by smooth in¬ 
variant curves for j3? Similarly, if ( S 2 ,g ) is a metric on S 2 which is isospectral to a surface 
of revolution, is S*S 2 foliated by smooth Lagrangian tori invariant under the geodesic flow? 
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